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Abstract 

It is shown that large classes of nonlinear systems of PDEs, with possi- 
bly associated initial and/or boundary value problems, can be solved 
by the method of order completion. The solutions obtained can be 
assimilated with Hausdorff continuous functions. The usual Navier- 
Stokes equations, as well as their various modifications aiming at a 
realistic modelling, are included as particular cases. The same holds 
for the critically important constitutive relations in various branches 
of Continuum Mechanics. The solution method does not involve func- 
tional analysis, nor various Sobolev or other spaces of distributions 
or generalized functions. The general and type independent existence 
and regularity results regarding solutions presented here have recently 
been introduced in the literature. 

"... provided also if need be that the notion of a solution 
shall be suitably extended ..." 

cited from Hilbert's 20th Problem 

1. Main ideas of the order completion solution method 

The solution method is divided in two parts. The proof of the existence 
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of solutions follows the method of order completion introduced and 
first developed in Oberguggenberger & Rosinger. The proof of the 
regularity of solutions is a consequence of recent results obtained in 
Anguelov [1] , regarding the structure of the Dedekind order completion 
of spaces of continuous functions C(X), where X is a rather arbitrary 
topological space. 

For simplicity of presentation, we shall consider single nonlinear PDEs. 
The extension to systems of such nonlinear PDEs and associated initial 
and/or boundary value problems can - rather surprisingly - be done 
easily, as seen in Oberguggenberger & Rosinger, this is being one of the 
major advantages of the order completion method. Let us therefore 
consider nonlinear PDEs of the general form 

(1.1) F(x,U(x), . . . ,D p x U(x), . . . ) = f(x), x eVl C R™ 

with p G N n , \p\ < m, where the domains Q can be any open, not nec- 
essarily bounded subsets of R™, while the orders m G N of the PDEs 
are arbitrary given, and the unknown functions, that is, the solutions 
one looks for are U : f2 — > R. 

The unprecedented generality of these nonlinear PDEs comes, above 
all, from the class of functions F which define the left hand terms, and 
which are only assumed to be jointly continuous in all of their argu- 
ments. The right hand terms / are also required to be continuous. 
As seen, however, both F and / can have certain discontinuities as 
well, Oberguggenberger & Rosinger. 

Regardless of the above generality of the nonlinear systems of PDEs 
considered, one can find for them solutions U defined on the whole of 
the respective domains Q. These solutions U have the blanket, type 
independent, or universal regularity property that they can be assimi- 
lated with Hausdorff continuous functions. 

It follows in this way that, when solving systems of nonlinear PDEs of 
the generality of those in (1.1), one can dispense with the various cus- 
tomary spaces of distributions, hyperfunctions, generalized functions, 
Sobolev spaces, and so on. Instead one can stay within the realms of 
usual functions, more precisely, of interval valued functions, see the 
Appendix for a short presentation of essentials on Hausdorff contin- 
uous functions. Also, when proving the existence and the mentioned 
type of regularity of such solutions one can dispense with methods of 
Functional Analysis. However, functional analytic methods can pos- 
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sibly be used in order to obtain further regularity or other desirable 
properties of such solutions. 

Let us now associate with each nonlinear PDE in (1.1) the correspond- 
ing nonlinear partial differential operator defined by the left hand side, 
namely 

(1.2) T(x, D)U(x) = F(x,U(x), . . . ,DPU(x), . . . ), x G Q 

Two facts about the nonlinear PDEs in (1.1) and the corresponding 
nonlinear partial differential operators T(x, D) in (1.2) are important 
and immediate : 

• The operators T(x, D) can naturally be seen as acting in the 
classical context, namely 

(1.3) T(x,D) : C m (n)3U i — ► T(x,D)U G C°(Q) 
while, unfortunately on the other hand : 

• The mappings in this natural classical context (1.3) are typically 
not surjective even in the case of linear T(x,D), and they are 
even less so in the general nonlinear case of (1.1), (1-2). 

In other words, linear or nonlinear PDEs in (1.1) typically cannot be 
expected to have classical solutions U G C m (Q), for arbitrary contin- 
uous right hand terms / G C°(fi), as illustrated by a variety of well 
known examples, some of them rather simple ones, see Oberguggen- 
berger & Rosinger [chap. 6]. 

Furthermore, it can often happen that nonclassical solutions do have 
a major applicative interest, thus they have to be sought out beyond 
the confines of the classical framework in (1.3). 

This is, therefore, how we are led to the necessity to consider gen- 
eralized solutions U for PDEs like those in (1.1), that is, solutions 
U ^ C m (Q), which therefore are no longer classical. This means that 
the natural classical mappings (1.3) must in certain suitable ways be 
extended to commutative diagrams 
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C m (Q) 



T(x,D) 



(1.4) 



C 



c 



X 



Y 



T 



with the generalized solutions now being found as 

(1.5) uex\c m (n) 

instead of the classical ones U G C m (Q) which may easily fail to exist. 
A further important point is that one expects to reestablish certain 
kind of surjectivity type properties typically missing in (1.3), at least 
such as for instance 

(1.6) C°{Q) C f(X) 

As it turns out, when constructing extensions of (1.3) given by com- 
mutative diagrams (1.4), we shall be interested in the following some- 
what larger spaces of piecewise smooth functions. For any integer 
< I < oo, we define 



;i-7) cm 



u : Q 



3 T C f2 closed, nowhere dense 
u e C l (n\T) 



and as an immediate strengthening of (1.3), we obviously obtain 



1.8) T(x,D)CZm C C° nd (n) 



The solution of the nonlinear PDEs in (1.1) through the order com- 
pletion method will come from the construction of specific instances 
of the commutative diagrams (1.4), given by, see (2.18), (2.27) 
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T(x,D) 



T 



(1.9) 



M°(n) 



bijective 



where, as elaborated later, the operation ( )# means the Dedekind 
order completion, according to MacNeille, of the respective spaces, as 
well as the extension to such completions of the respective mappings, 
see Oberguggenberger & Rosinger [Appendix]. It follows that in terms 
of (1.4), we have 



thus we shall obtain for the nonlinear PDEs in (1.1) generalized solu- 
tions 



Furthermore, instead of the surjectivity condition (1.6), we shall at 
least have the following stronger one 



So far about the main ideas related to the existence of solutions of 
general nonlinear PDEs of the form (1.1). 

As for the regularity of such solutions, we recall that, as shown in 
Oberguggenberger & Rosinger, one has the inclusions 



where Mes (Q) denotes the set of Lebesgue measurable functions on Q. 
In this way, in view of (1.9), (1.10), one can assimilate the generalized 
solutions U of the nonlinear PDEs in (1.1) with usual measurable 
functions in Mes (Q). 



X = M%{Q)#, Y = M°(n)#, T = T# 



(1.10) UeM^iSl)* 



(1.11) C(fi) C T#(MPM#) 



(1.12) M°{Q)* C Mes{Q) 
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Recently, however, based on results in Anguelov [1] , it was shown that 
instead of (1.12), one has the much stronger property 

(1.13) M°{Q)# C H(fi) 

where H (Q) denotes the set of Hausdorff continuous functions on fl. 
Consequently, now one can significantly improve on the earlier regu- 
larity result, as one can assimilate the generalized solutions U of the 
nonlinear PDEs in (1.1) with usual functions in H(Q). 



2. The construction of diagram (1.9) 

Since we solve PDEs through order completion, let us see how near we 
can come to satisfying the equality in (1.1), when using inequalities. 
For that purpose, it is useful to consider for each x G £1 the following 
set of real numbers 

(2.1) R x = { F{x,£ , ...,&,...) | & e R,forp e N n , \p\<m} 

Clearly, for x E Q fixed, R x is the range in R of F(x, . . . ), and 
since F is jointly continuous in all its arguments, it follows that R x is 
a nonvoid interval which is bounded, half bounded, or is the whole of 
1BL This latter case, which can happen often with nonlinear PDEs in 

(1.1) , will be easier to deal with, see (2.3) next. Clearly, in the case of 
non-degenerate linear PDEs in (1.1), this latter case always happens. 
Given now x G Q, it is obvious that a necessary condition for the 
existence of a classical smooth solution U G C m of (1.1) in a neigh- 
bourhood of x is the condition 

(2.2) f(x) G R x 

Consequently, for the time being, we shall make the assumption that 
the right hand term functions / in the nonlinear PDEs in (1.1) satisfy 
the somewhat stronger version of condition (2.2) given by 

(2.3) f(x) G int R x , for x G Q 
Clearly, whenever we have 

(2.4) R x — R, for x G tt 
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then (2.3) is satisfied. And as mentioned, this is the case with all 
nontrivial linear PDEs, as well as with most of the nonlinear PDEs of 
practical interest. 

And now the basic and rather simple local approximation result on 
how near we can satisfy the equality in (1.1), when using inequalities. 



Proposition 2.1. 

Given / G C°(ft), then 

V x G ft, e > : 

?2 z\ 3 5 > 0, P polynomial in x G R n : 

V x G ft, ||rr — rr || < 5 : 
f(x)-e < T(x,D)P(x) < f(x) 

Proof 

Given x G ft, then for e > small enough, condition (2.3) yields 
£ p G R, with p G N n , \p\ < m, such that 

(2.6) F(x ,£ , . . ..C, ) f(x Q )-e/2 

Let us take P a polynomial in x G R n , which satisfies the conditions 

DP x P{x ) = C P , p G W\ \p\ < m 

In this case from (2.6) we clearly obtain the relation 

(2.7) T(x ,D)P(x ) = f(x )-e/2 

and since both T(x, D)P(x) and f(x) are continuous in x G Q, the 
local inequality property (2.5) follows easily from (2.7). 

The global approximation version of the inequality property in (2.5) is 
given in 

Proposition 2.2. 

If / G C°{Q), then 

V e > : 

(2.8) 3 T e C tt closed, nowhere dense, U £ G C m {Vt \ Y e ) : 

f-e < T(x,D)U € < f on ft \ T e 



Proof 

Let us take a covering of Q of the form 



(2.9) n 




V 



where K v are compact n-dimensional intervals in M n , namely, K v = [a v , b u ], 
with a v = (a Vil , . . .,a Vin ), b u = (b Vjl , . . .,b ViTl ) G M n . We also assume, 
see Forster, that the covering (2.9) is locally finite, that is 

V xeVl : 

(2.10) 3 V x C Q neighbourhood of x : 

{ i/ E N \ K v C\V X ^ (j> } is a finite set of indices 

and furthermore 

(2.11) the interiors of K u , with v e N, are pairwise disjoint 

Let us now take e > arbitrary but fixed. Further, we take v e N. We 
shall apply Proposition 2.1 to each xq G K v . Then we obtain 5 X0 > 
and a polynomial P XQ such that 

/(x)-e < T(x,D)P X0 (x) < f(x), xen, \\x-x \\ < 5 Xo 

But K v is compact, therefore 

3 5 > : 

V x G K v : 

(2.12) 3 P X0 polynomial in x G lR n : 



Now we shall subdivide K u , which was assumed to be a compact n- 
dimensional interval, into n-dimensional subintervals I±, . . . so 
that the diameter of each of them is less or equal 5. 
Let us denote by J a generic such n-dimensional subinterval in any of 
the K u , when v G N. If aj G J is the center of any such n-dimensional 
subinterval then (2.12) gives for x G int J 




f(x)-e < T(x,D)P X0 (x) < f(x) 



f(x)-e < T(x,D)P aj (x) < f(x) 



Let us now take 



(2.13) r e = D, \ Uj int J 
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that is, with the union ranging over all such n-dimensional subintervals 
J. In view of ??? If we define U e E C m {VL \ Y e ) by 

U e = P aj on p| int J 

then the proof is completed. 

Remark 2.1 

1) It is easy to see that the inequalities in (2.5) and (2.8) can be 
replaced with the following ones, respectively 

(2.14) f(x) < T(x,D)P(x) < f(x)+e 

(2.15) / < T(x, D)U e < f + e 

as their proofs follow after the corresponding obvious minor changes in 
the proofs of the above two propositions. And these four inequalities 
are sharper than would respectively be the inequalities 

f(x)-e < T(x,D)P(x) < f(x) + e, f-e < T(x,D)U e < f + e 

As we shall see not much later, we do need the sharper inequalities. 
Indeed, the order completion method which we shall employ is based 
on MacNeille's construction, therefore, it uses Dedekind cuts. And 
such cuts do need the above sharper inequalities. 

2) In Proposition 2.2., as well as in its version corresponding to the 
above inequality (2.15), we have in addition the property 

(2.16) mes (T e ) = 

where mes denotes the usual Lebesgue measure. Indeed, according to 
(2.10), (2.11) and (2.13), T e is a countable union of rectangular grids, 
each generated by a finite number of hyperplanes. 
Here it should be noted that the presence of the closed, nowhere dense 
singularity sets T € in the global inequalities (2.8) and (2.15) will prove 
not to be a hindrance. And in fact, it will lead to the classes of piece- 
wise smooth functions in (1.7) which prove to be convenient. 
The presence of such closed, nowhere dense singularity sets is rather 
deep rooted, as it is connected with such facts as the flabbiness of re- 
lated sheaves of functions, see Oberguggenberger & Rosinger [chapter 
7], or the global version of the classical Cauchy-Kovalevskaia theorem 
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on analytic nonlinear PDEs, see Oberguggenberger & Rosinger and 
the literature cited there. 

3) As seen from the proof of Proposition 2.2., the functions U e can in 
fact be chosen as piecewise polynomials in x G M™. 

Let us now note that there is an obvious ambiguity with the piecewise 
smooth functions in C l nd (Q) in (1.7). Indeed, given any such function 
u G C l nd (Q), the corresponding closed, nowhere dense set T cannot 
be defined uniquely. Therefore, it is convenient to factor out this 
ambiguity, and this can be done easily as follows. Since C® d (Q) is 
the largest of these spaces of functions, we shall do for this space 
the mentioned factoring out, by defining on it the equivalence relation 
u ~ v for any two elements u, v G C° d (fi), as given by the condition 

3 TcO closed, nowhere dense : 

(2.17) *) u, v g c°(n\r) 

* *) u = v on n\r 

It is easy to see that m defined above is indeed an equivalence 
relation, since the union of a finite number of closed and nowhere dense 
subsets is again closed and nowhere dense. Now we can eliminate the 
mentioned ambiguity by going to the quotient space 

(2.18) M°(Q) = M/ « 

and in view of (2.8), (2.15), we define for any two elements G, H G 
M°(£l) the partial order G < H , by 

3 g E G, h G H, rcfl closed, nowhere dense : 

(2.19) *) g, he C°(n\T) 

* *) g < h on Q, \ T 

Let us now denote by 

(2.20) (M°(n)#,<) 

the Dedekind order completion due to MacNeille of the partially or- 
dered space (M°(ft), <) which was defined in (2.18), (2.19). Then this 
space A4°(Q) # in (2.20) is order complete, and we have the order iso- 
morphical embedding, see Oberguggenberger & Rosinger [Appendix] 

(2.21) M°(Q) 3 G i — > <G]e M°(n)* 
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which also preserves the infima and suprema. Moreover, in view of the 

MacNeille construction, we can further extend the embedding (2.21) 

as follows . . , 

o. 1. e. id 

M°(Q) >- M°(n)# >- V(M°(tt)) 

(2 ' 22) G . <G] *<G] 

In order to obtain the full situation with respect to the range of the 
nonlinear partial differential operators T(x, D), we note that C°(Q) C 
C® d (Q), and we have the order isomorphical embedding 

(2.23) C°{n) 3g \ — > Ge M°{n) 

where G is the ps equivalence class of g. Furthermore, the partial 
order < on induces on C°(Q) through this embedding the 

usual point-wise order relation of functions, namely, g < h, if and only 
if g(x) < h(x), for x G Q. 

Let us now recall that our main interest is the construction of the 
commutative diagrams (1.9). In this regard, having constructed in 
(2.18) - (2.20), respectively, the spaces and M°(tt)* and their 

partial orders, the next step is to construct the partially ordered spaces 
M.™{£1) and M^(VL)* . For that purpose we start with (1.11), namely 

(2.24) T(x,D):CZm — > C° nd (Q) 

As mentioned, we shall solve the nonlinear PDEs in (1.1) by extending 
through order completion this mapping in (2.24), and we do so by 
constructing the commutative diagrams in (1.9). And at this stage we 
are now in the position to start doing so step by step. Let us note 
first that if u e C m (Q \ T), where r C Q is any given closed, nowhere 
dense subset, then we also have 

(2.25) T(x,D)ueC°(n\T) 

This means that the singularity subsets T do not increase by the appli- 
cation of the nonlinear partial differential operators T(x, D). However, 
as before with C® d (Q), the ambiguity about associating such singular- 
ity sets to functions in C^(Q) remains. Therefore, in view of (2.25), 
we shall define the equivalence relation u «t v, for any two functions 
u, v e C^(fi), by the condition 

(2.26) T(x,D)u w T(x, D)v 
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which uses the equivalence relation given in (2.17), and in 

addition, it also depends on the nonlinear partial differential operator 
T(x,D). In fact, this equivalence relation on C™ d (Q) is what 

is called the pull-back through the mapping T(x, D) in (2.24) of the 
equivalence relation on C° d (fi). 
Let us now define the quotient space 

(2.27) A<5?(n) = cs(n)/ 

in which case the mapping (2.24) generates canonically the injective 
mapping 

(2.28) T:M%{n) — ► A4°(fi) 

defined by T(U) = G, where G is the unique rs equivalence class in 
.M°(f2) of any of the T(x, -D)w, where u belongs to the m T equivalence 
class U in M^(Q). 

At last, we can define the partial order on M.™{TL) as the 

pull-back through the mapping T in (2.28) of the partial order < 
in (2.19) on M°(n), that is, for U, V e M%(tt), we have U < T V, if 
and only if 

(2.29) T U < TV 

In this way we obtain the partially ordered set (.M™(f2), <t) giving 
the desired order structure on the domain A4Jp(Q) of T, which is 
the mapping in (2.28) that corresponds now to our nonlinear partial 
differential operator T(x, D) in (1.2), (1.3), or more precisely, in (2.24). 
It is obvious in view of (2.29) that the injective mapping T in (2.28) 
is also an order isomorphical embedding. 

So far, we have in this way obtained the top commutative rectangle 
in (1.9). 

Applying now to (A4™(Q), < T ) the Dedekind order completion of Mac- 
Neille, we obtain 

(2.30) (MJp(n)*,< T ) 

which is order complete, and in addition, similar with (2.21), we also 
have the order isomorphical embedding 

(2.31) M^(tt) 3 U i — ► <U]e M^(tt)* 
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which preserves the infima and suprema. Also, similar with (2.22), we 
have 



o. i. e. id 

A4!p(fi) M^{9)# V(M%(n)) 

(2-32) 

V ; U *■ <U\ * <U\ 



And now all that remains is to define T* in (1.9). In view of (2.32), 
however, this can be done in a standard manner following from the 
MacNeille order completion, see Oberguggenberger & Rosinger [Ap- 
pendix]. Consequently, one obtains the order isomorphical embedding 

(2.33) T* : M%(n))# — > M°(Q)* 

which also preserves the infima and the suprema. In more detail, we 
have the following commutative diagram 



T 

M%(Q.) 3 U ► T(U) E M°(n) 

(2.34) 



M^{Q)*3 <U\ ► T#{<U\) = 

T* = < T(U) ] E M°(n) 

In this way we have indeed obtained the whole of the commutative 
diagram (1.9), which we shall present now in the form seen next. Here 
"sur" and "inj" mean mappings which are surjective and injective, 
respectively, while as before, "o. i. e." means order isomorphical em- 
bedding, and "o. i." stands for order isomorphism. The dotted arrows 

"< " mean the "pull-back" through which the respective 

structures were defined 
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T(x,D) 



sur 



i T < 



sur 



(2.35) 



•'- / < 



o. 1. e. 



< 



mj 



o. i., see (3.1) below 
M^{VL)# ► M°(Q)* 



3. General existence result 



One of the typical main existence results concerning the solutions of 
the nonlinear PDEs in (1.1) is presented in the following theorem, see 
Oberguggenberger & Rosinger [38-64] for a proof 



Theorem 3.1. 



(3.1) T# 



(Q) 



# 



M°(Q) 



# 



This means that, given the nonlinear PDEs in (1.1), for every right 
hand term / G A4°(Q)^, there exists a generalized solution U G 
M.™ (f2) # , satisfying the relation T#U = f, according to the extension 
in (1.9). 

As seen in Oberguggenberger & Rosinger [74-93], the space A4°(Q)# 
in which the right hand terms / of the nonlinear PDEs in (1.1) can 
range - and which now are solved by Theorem 3.1. - contains a large 
amount of discontinuous function on Q. Certainly, in view of (2.18), 
_M°(fi)# contains all the piecewise discontinuous functions in C® d (Q). 
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What is particularly important to note is that, in view of (3.1), a 
variety of linear and nonlinear PDEs can be solved, in spite of the 
fact that the respective PDEs are known not to have solutions in 
distributions. Among them is the celebrated 1957 Hans Lewy example, 
see Obrguggenberger & Rosinger [chap. 6, 8]. 

In this regard, it was in Oberguggenberger & Rosinger that this Hans 
Lewy example of a PDE not solvable in distribution was nevertheless 
solved for the first time through the method of order completion. 

The coherence between the solutions obtained in (3.1) and the usual 
classical solutions, whenever the nonlinear PDEs in (1.1) may have the 
latter, follows easily from the commutative diagram (2.35). In other 
words, whenever the nonlinear PDEs in (1.1) happen to have classical 
solutions U E C m (f2), then they are also generalized solution in the 
sense of (3.1). 

Finally, it is important to note that the above existence result in (3.1) 
can easily be extended to systems of nonlinear PDEs of the general 
form in (1.1), see Oberguggenberger & Rosinger [chap. 8-11]. 



4. Initial and/or boundary value problems 
and constitutive relations 

One of the significant advantages of the order completion method in 
solving systems of nonlinear PDEs of the general form in (1.1) comes 
from the ease initial and/or boundary value problems associated with 
such equations can be solved. This is in strong contradistinction with 
the variety of functional analytic methods of solution where consid- 
erable difficulties arise related to the need to restrict distributions or 
generalized functions to lower dimensional manifolds. Indeed, such 
operations of restriction are typically ill-defined. 

On the other hand, when using the method of order completion, the 
issue of satisfying the initial and/or boundary values can be decou- 
pled from the issue of the existence of solutions. Indeed, satisfying the 
initial and/or boundary values can be dealt with first and separately 
from the issue of proving the existence of solutions. 
Details in this regard can be found in Oberguggenberger & Rosinger 
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[chap. 8, 11]. And the reason behind that rather surprising ease the 
order completion method exhibits when dealing with initial and/or 
boundary value problems comes from a fact seen next, in section 5. 

In Fluid Dynamics, and in general, Continuum Mechanics, a critical 
role is played by constitutive equations, see Rajagopal & Wineman, 
Rajagopal, or Rajagopal & Srinivasa. 

The usual functional analytic methods can - if at all - deal with such 
constitutive equations in no less a difficult manner than they can do 
with initial and boundary value conditions. A regrettable consequence 
of these considerable difficulties is the failure so far of functional an- 
alytic methods to approach in any significant, let alone, systematic 
manner, the issue of such critically important constitutive equations. 
Here again, the order completion method proves its advantage by be- 
ing able to deal as well with constitutive equations. This is simply 
a consequence of the fact that, as mentioned next in section 5, the 
order completion method can solve equations which are far more gen- 
eral than the linear or nonlinear systems of PDEs, or the constitutive 
equations. 



5. An abstract existence result 

A better understanding of the power underlying the order comple- 
tion method in solving very large classes of equations, classes far be- 
yond any nonlinear systems of PDEs, can be obtained from the follow- 
ing rather abstract existence result, see Oberguggenberger & Rosinger 
[chap. 9]. 

Let X be any set, and let (Y, <) be any partially ordered set which 
has no minimum or maximum. Further, let 

(5.1) T : X — > Y 

be any given mapping. The problem we consider is to find a solution 
A G X for the equation 

(5.2) T(A) = F 
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for any given FgF. The answer is obtained as follows. Similar with 
the construction of the commutative diagrams (1.9) and (2.35), one 
can construct commutative diagrams 



T 



X 



Y 



T 



(5.3) 



Xj> 



Y 




bijective 



Y* 



And then the following result on the existence of solutions holds 
Theorem 5.1. 

For any given F G Y#, the equation 



The significant generality of the above existence result allows, among 
others, the separation mentioned in section 4, between first satisfying 
the initial and/or boundary value conditions, and then, second, prov- 
ing the existence of solutions in the case of general nonlinear systems of 
PDEs of the form in (1.1). Indeed, by first imposing the initial and/or 
boundary values, one is in fact defining the set X in (5.1). And as 
seen above, that can be done without any restrictions. Subsequently, 
condition (5.5) is both necessary and sufficient for the existence of a 
generalized solution A e X* . 
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6. The Hausdorff continuity of solutions 

The major novelty in this paper is about the fact that the solutions 
U G .M™(f2)* of systems of nonlinear PDEs of type (1.1), obtained 
according to the procedure in Theorem 3.1., and of its generalizations 
can now be assimilated with Hausdorff continuous functions in H(Q). 
In fact, as seen in (A. 12) in the Appendix, the mentioned solutions 
can be assimilated with nearly finite Hausdorff continuous functions. 

7. Final comments 

A further advantage of the order completion method is that one is not 
limited to consider in (1.9) and (2.35) only the pull-back partial order 
<t generated by the partial differential operators T(x, D) on 
Indeed, as seen in Oberguggenberger & Rosinger [chap. 13], a large 
variety of other partial orders on M^(Vt) can still secure existence 
theorems similar to Theorem 3.1. 

As for the use of pull-back orders, it is important to note that there 
exists a certain analogy with functional analytic methods for solving 
PDEs. Indeed, in such methods, the topologies considered on the do- 
mains of the partial differential operators T(x, D) are but pull-backs 
through these operators of suitable topologies on their ranges. 
Details about such facts, and in general, about certain similarities be- 
tween the order completion method and the usual functional analytic 
ones in solving PDEs can be found in Oberguggenberger & Rosinger 
[chap. 12]. 

The results in this paper invite a comparison with the customary per- 
ception regarding the solution of linear or nonlinear PDEs. Typical 
for that perception are the following to recent citations. 

The 2004 edition of the Springer Universitext book " Lectures on PDEs" 
by V I Arnold, starts on page 1 with the statement : 

"In contrast to ordinary differential equations, there is no 
unified theory of partial differential equations. Some equa- 
tions have their own theories, while others have no theory 
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at all. The reason for this complexity is a more compli- 
cated geometry ..." (italics added) 

Similarly, the 1998 edition of the book "Partial Differential Equations" 
by L C Evans, starts his Examples on page 3 with the statement : 

"There is no general theory known concerning the solv- 
ability of all partial differential equations. Such a theory 
is extremely unlikely to exist, given the rich variety of phys- 
ical, geometric, and probabilistic phenomena which can be 
modelled by PDE. Instead, research focuses on various par- 
ticular partial differential equations ..." (italics added) 

Appendix : Definition and Properties of 

Hausdorff-Continuous Functions 

The Hausdorff continuous functions are not unlike the usual real val- 
ued continuous functions. For instance, they assume real values on a 
dense subset of their domain of definition and are completely deter- 
mined by the values on this subset. However, these functions may also 
assume interval values on a certain subset of their domain of defini- 
tion. Hence the concept of Hausdorff continuity is formulated within 
the realm of interval valued functions. We shall deal in this Appendix 
with functions whose values can be not only usual real numbers but 
also extended real numbers, that is, elements in R = R U {— oo, +00}. 
Moreover, as mentioned it proves to be convenient to allow the values 
of the functions to be not only numbers in R, but also closed intervals 
of such numbers, namely, [a, b] C R, with a, b G R, a < b. 

Towards the end of the 19-th century, Baire brought in the concepts 
of lower and upper semi- continuous functions, when dealing with non- 
smooth real valued functions. And in effect, he associated with each 
real valued function /, two other real, or extended real valued functions 
/(/) and S(f), with /(/) < / < S(f), which proved to be particularly 
helpful, see (A. 6), (A. 7). However, following the prevailing mentality 
at the time, each of these three functions was considered separately 
and as being a single valued function. 
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As it turns out on the other hand, by considering interval valued func- 
tions, such as for instance F(f) = [1(f), S(f)], one can significantly 
improve on the understanding and handling of non- continuous func- 
tions. 

The study of interval valued functions can, among others, show that 
the particular case of functions which have values given by one single 
number is appropriate for continuous functions only. On the other 
hand, non-continuous functions are much better described by suitably 
associated interval valued functions. 

Indeed, in the case of functions / which are not continuous, a much 
better description can be obtained by considering them given by a 
pair of usual point valued functions, namely / = [/,/], thus leading 
to interval valued functions, according to f(x) = [f(x),f(x)] C R, 
for x G Q. And then, a natural class which replaces, and also ex- 
tends, the usual point valued continuous functions is that of Hausdorff- 
continuous interval valued functions, see below Definition Al. The 
distinctive and essential feature of these Hausdorff-continuous func- 
tions / = [ /, / ] is a condition of minimality with respect to the gap 
between / and /, with the further requirement that / be lower semi- 
continuous, and / be upper semi-continuous. 

The interest in more recent times in interval valued functions comes 
from a number of branches of mathematics, such as approximation 
theory, Sendov, and numerical analysis, Kraemer. 

Most of the results on interval valued functions presented in this Ap- 
pendix have, however, been developed by R Anguelov. This was done 
in view of the usefulness of such functions in several branches of math- 
ematics, see Anguelov [1,2], Anguelov & Markov, Anguelov et.al [1,2], 
Anguelov & Minani, Anguelov & Rosinger [1-3]. Here, owing to re- 
striction of space, we shall only present a minimal amount of them, 
needed in order to support section 6 above. For the full details in this 
regard, including proofs, see Appendix 2 in Anguelov & Rosinger [3]. 

As mentioned, the class of interval valued functions of special inter- 
est here is that of Hausdorff-continuous, or in short, H-continuous 
functions. As it turns out they enjoy the minimality property (A. 8) 
with respect to their graph completion, see (A. 7), and that allows an 
effective interplay between Analysis and Topology, with the latter in- 
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volving both the domain and the range of the functions dealt with. 
Let 



(A.l) IR = { [a, a] | a, a G R = R U {— oo, +00}, a < a} 

be the set of all finite or infinite closed intervals. The functions which 
we consider can be defined on arbitrary topological spaces Q. For the 
purposes of the nonlinear PDEs studied in this paper, however, it will 
be sufficient to assume that Q C R n are arbitrary open subsets. Let 
us now consider the set of interval valued functions 

(A.2) A(n) = {f-.n — > iR } 

By identifying the point a G R with the degenerate interval [a, a] G IR, 
we consider R as a subset of IR. In this way A(f2) will contain the set 
of functions with extended real values, namely 

(A.3) A{Q) = { f : Q — > R } C A(Q) 

We define a partial order < on IR by 

(A.4) [a, a] < [b, b] <^ a < b, a < b 

Now on A(f2) we define the partial order induced by (A.4) in the usual 
point- wise way, namely, for f,g<E A(fi), we have 

(A.5) f<g ^ f(x)<g(x), x£ll 

Clearly, when restricted to A(Q), the above partial order on A(f2) 
reduces to the usual one among point valued functions. 

Let / G A(f2). For every x G f2, the value of / is an interval, 
namely, f(x) = [ f(x), J(x) ], with f_(x), J(x) G R, f_(x) < 
f(x). Hence, every function / G A(Q) can be written in the form 
/ = [/,/], with /, /_G ^(fl), / < / < /, and 

/ e A(Q) ^ 1 = f = f. 

In the particular case of functions in A(fl), that is, with extended real, 
but point, and not nondegenerate interval values, a number of basic 
results were obtained already in Baire, see also Nicolescu for a more 
recent detailed presentation. The rest of the more general results con- 
cerning functions in A(f2), that is, with values finite or infinite closed 
intervals, were developed for the first time in the above cited works 
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of Anguelov . The few such earlier results were obtained in Sendov, 
where the particular instance of ft C R was dealt with. 

For x G ft, we denote by Bs (x) the open ball of radius 5 centered at x. 
Let us consider any dense subset .D C ft, and associate with it the pair 
of mappings I(D,Q,.), S(D,Q,.) : A(ft) — > A(fl), called /operand 
upper Baire operators, respectively, where for every function / G A(ft) 
and x G ft, we define 

I(D,Q,f)(x) = sup 5>0 inf {^G/(y) | y G B$(x) D D } 
{ - 0) S(D,Q,f)(x) = mi s>0 sup {zef(y) \ y G B$(x) (1 D } 

In Baire, these two operators were considered and studied in the par- 
ticular case of functions / G A(£l) and when D — Q, see also Nicolescu. 
In view of the main interest here in interval valued functions / G A(Q), 
it is useful to consider as well the following third mapping, namely, 
F(D, tt, f) : A(Q.) -> A(Q), defined for / G A(Q) by 

(A.7) F(D,nj)(x) = [i(D,nj)(x), s(D,nj)(x)\, xen, 

and called the graph completion operator. In case D = Q, we use the 
simpler notations /(ft, ft, /) = /(/), 5(ft, ft, /) = S(f) and F(ft, ft, /) = 

^(/)- 

The next definition was given in Sendov in the case of ft C R, however, 
it can obviously be extended to any topological space ft. 

Definition A.l 

A function / G A(ft) is called Haus dor -ff- continuous, or in short, H- 
continuous, if and only if for every function g G A(ft), we have satisfied 
the following minimality condition on / 

(A.8) g(x) C f(x), xeQ =}► F{g) = f 

We shall denote by H(ft) the set of all Hausdorff-continuous interval 
valued functions on ft. 

Definition A2 

A function / G H(ft) is called nearly finite, if and only if there exists 
an open and dense subset D C ft, such that 

(A. 9) f(x) G IK is a finite interval for x G D 
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We denote by M n f(il) the set of nearly finite H-continuous functions 

Regarding the regularity properties of solutions of general nonlinear 
systems of PDEs of the form in (1.1), a crucial role is played by the 
following mapping 

(A.10) F :C° d (n)3u i— > F(fl\r ) fl ) «)6H l/ (fl) 

where we recall that, according to (1.7), for every u G C° d (f2), there 
exists a closed, nowhere dense subset rcfi, such that u G C°(f2 \ T), 
hence in view of (A. 7), F(Q \ T,Q,u) is well defined. The fact that 
such a T need not be unique, does not affect the above definition, see 
Anguelov & Rosinger [3]. 

The following theorem shows that the images of two functions in 
C n d(Q) under the mapping F in (A. 10) are the same, if and only if 
these functions are equivalent with respect to the equivalence relation 
(2.17). 

Theorem A.l 

Let u, v G C nd (Q). Then F (u) = F (v) -<=^ u m v 

In view of (2.18), (A. 10) and the above theorem now we can define a 
mapping 

(A.ll) F : M°(n) — >m nf (n) 

in the following way. Let u G U G M°(£i), then F (U) = F(u). 
It is easy to see that the definition of F Q (U) does not depend on the 
particular representative u G U of the equivalence class U. 

Theorem A. 2 

The mapping F : Ai°(Q) — > M n f(Q) defined in (A.ll) is an or- 
der isomorphic embedding with respect to the order relation (2.19) on 
and the order relation induced by (A. 5) on H n /(f2). Namely, 
for any U, V G M°(Q), we have 

U < V ^ F (U) < F Q (V) 
Finally, we also have 
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Theorem A. 3 

The set H n /(f2) is Dedekind order complete with respect to the partial 
order induced on it by (A. 5). 

Let g G M. n f(Q). Then there exists a subset Q C JH°(Q) such that 
g = su P F (£) = sup{F (G) | G E Q} 

This theorem shows that H n /(f2) is the smallest Dedekind order com- 
plete subset of HI(f2) which contains the image of Ai°(£l) under the 
order isomorphical embedding Fq. Hence it is order isomorphic to the 
Dedekind order completion Ai°(fl)^ of In this way we obtain 

the commutative diagram, where F # denotes the order isomorphism 
from M°(n)* to H n/ (fi), namely 

C° n M —5- M nf (Q) 

graph completion ■ 

order isomorphical embedding ■ 

M°(tt)* Hn/(n) 

order isomorphism 

Now we can bring together the above diagram with the one in (2.35) 
and obtain the two successive order isomorphisms 

T# F* 
(A12) M^(Sl))* *- M°(n)* 2 ^ H n/ (Q) 

It follows therefore that the set in which the solutions of 

the general nonlinear systems of PDEs of the form in (1.1) are found 
is mapped by the bijection Ff o T# onto the set M. n f(Q) of all nearly 
finite Hausdorff continuous functions. Since both these mappings are 
order isomorphisms, the set A4^(Q))^ is order isomorphic with the 
set H n/ (Q). 

Hence, the solutions of the general nonlinear systems of PDEs of 
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the form in (1.1), which are obtained through the order completion 
method, can always be assimilated with nearly finite Hausdorff con- 
tinuous functions. 

This is the argument supporting section 6 above. 

For the sake of further clarification, let us turn now in the remaining 
part of this Appendix to some of the issues concerning the disconti- 
nuities of Hausdorff-continuous functions. Arbitrary interval valued 
functions / = [ /, / ] e A(Q) can exhibit a variety of types of discon- 
tinuities, and certainly not less so, than usual point valued functions 
/ E A(D.) do. 

Hausdorff-continuous functions, although generalize usual point val- 
ued continuous functions, enjoy nevertheless a number of nontrivial 
continuity related properties. On the other hand, Hausdorff-continuous 
functions can have quite large sets of discontinuities. This shows that 
they do indeed form a larger class than the usual continuous functions, 
even if they still have important similar properties. 
The fact that large enough sets of discontinuities can be present with 
Hausdorff-continuous functions allows for their use - as seen in this 
paper - in obtaining the existence of nonclassical solutions for large 
classes of systems of nonlinear PDEs. 

Consequently, and as mentioned, Hausdorff-continuous functions - pre- 
cisely since they are not generalized functions - can be seen as setting 
aside to a certain extent the variety of distributional and other tradi- 
tional generalized solutions of linear and nonlinear PDEs which have 
been obtained by functional analytic methods, or by the methods of 
the nonlinear algebraic theory listed by the AMS Subject Classifica- 
tion 2000, under 46F30. 

Indeed, solving large classes of nonlinear PDEs through Hausdorff- 
continuous functions offers, among others, the following double ad- 
vantage : 

• one can bring in a significant simplification by avoiding the va- 
riety of usual functional analytic methods with their spaces of 
distributions or generalized functions, and 

• one can obtain universal regularity results for solutions of large 
classes of systems of nonlinear PDEs. 
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Theorems A. 4 and A. 5 below show important properties of Hausdorff- 
continuous functions related to their sets of discontinuities. For ev- 
ery interval valued function /=[/,/]€ A(VL), we denote by 
r(/) = { x G VL ; f(x) < f(x) }, which is the set of points 
x G VL where / assumes values f(x) = [ f(x), f(x) ] that are non- 
degenerate intervals, and not merely points. 

It follows that at points x G r(/), the interval valued function / = 
[ /, / ] cannot be continuous in the usual sense, since it is not a usual 
point valued function. 

In the particular case of Hausdorff- continuous functions, this fact can 
further be clarified. Namely, given any point x G VL, then 

, , x G r(/) •<=>- / and / not continuous at x •<=>- 

-<=>- / or / not continuous at x •<=>- f(x) < f(x) 

And now the basic result on the discontinuities of Hausdorff-continuous 
functions 

Theorem A. 4 

Given any H-continuous function / G A(VL). Then r(/) is of first 
Baire category in VL. 

Let us further specify the structure of the discontinuity set r(/). For 
e > 0, let us denote r e (/) = { x G VL \ J(x) - fix) > e }. Then 
clearly = (J e >o r e(/) = U„>i r 1/n (/). 

The next theorem gives a further insight into the structure of the 
discontinuity set r(/) of Hausdorff-continuous functions / G A(VL). 

Theorem A. 5 

If the function / G A(VL) is H-continuous, then for every e > 0, the set 
r e (/) is closed and nowhere dense in VL. 

An important similarity between usual continuous, and on the other 
hand, Hausdorff-continuous functions is that both of them are deter- 
mined uniquely if they are known on a dense subset of their domains of 
definition. This property comes in spite of the fact that, as seen in the 
previous section, Hausdorff-continuous functions can have discontinu- 
ities on sets of first Baire category, and such sets can have arbitrary 
large positive Lebesgue measure, see Oxtoby. Indeed, we have 
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Theorem A. 6 



Let f = [ f, f ], g = [ g, g ] e A(fi) be two H-continuous functions, 
and suppose given any dense subset D C f2. Then 



The real line R is Dedekind order complete, but not order complete, 
while the extended real line R is both Dedekind order complete and 
order complete. 

Let us recall that a partially ordered set which is order complete will 
also be Dedekind order complete, but as seen above, not necessarily 
the other way round as well. 

Typically, various spaces of real valued functions encountered in Anal- 
ysis are neither Dedekind order complete, nor order complete, when 
considered with the natural point-wise partial order relation. 
However, as we can see next, this situation changes when we deal with 
the set M(Q) of Hausdorff-continuous functions. 

Theorem A. 7 (Anguelov [2]) 

The set HI(f2) of Hausdorff-continuous functions is order complete 
when considered with the partial order in (A. 5). 

The Dedekind order completeness of the space H(f2) of Hausdorff- 
continuous functions is a nontrivial property, in view of the various 
connections between the usual continuous, and on the other hand 
Hausdorff-continuous functions. Indeed, the space C(O) of usual real 
valued continuous functions, which we have seen is strictly contained 
in H(Q), is well known not to be Dedekind order complete. On the 
other hand, once H(f2) proves to be Dedekind order complete, its or- 
der completeness follows easily from the fact that R is order complete. 
Indeed, the smallest and largest elements in H(f2) are respectively the 
functions Q 3 x i — ► — oo and Q 3 x i — > +oo. 

As is well known and shown by simple examples the spaces of real 
valued continuous functions C(Q) are not Dedekind order complete, 
thus, are not order complete either. 

Since these spaces are partially ordered in a natural way, one can apply 
to them the MacNeille version of Dedekind order completion method, 
see Oberguggenberger & Rosinger [Appendix], Luxemburg & Zaanen, 



/(*) 
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or Zaanen. 

This however being a general construction based on Dedekind type 
cuts, it leaves open the question of the nature of the elements which 
are added to these spaces of continuous functions by the respective 
Dedekind order completion process. 

A classical, 1950 result in this regard was obtained by Dilworth in the 
case of bounded and real valued continuous functions on Q for arbi- 
trary completely regular topological spaces Q. Namely, the respective 
Dedekind order completion is given by all the normal upper semi- 
continuous functions on Q. Regarding the Dedekind order completion 
of spaces C(Q) certain results were obtained in Mack & Johnson. 

As seen in Anguelov [1], the Dedekind order completion of spaces C(Q) 
of real valued continuous functions was for the first time effectively 
constructed for a large class of topological spaces Q. In this construc- 
tion Hausdorff-continuous functions and some of their subspaces play 
a crucial role. Consequently, both the problem of the completion of 
the spaces C(fl), as well as that of the structure of the elements which 
are added to these spaces C{Q) through the respective completion find 
a convenient solution through the use of interval valued functions. 
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